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Abstract 

By extending the concept of Euler-angle rotations to more than three di- 
mensions, we develop the systematics under rotations in higher-dimensional 
space for a novel set of hyperspherical harmonics. Applying this formalism, 
we determine all pairwise Coulomb interactions in a few-body system with- 
out recourse to multipole expansions. Our approach combines the advantages 
of relative coordinates with those of the hyperspherical description. In the 
present method, each Coulomb matrix element reduces to the "1/r" form fa- 
miliar from the two-body problem. Consequently, our calculation accounts 
for all the cusps in the wave function whenever an inter-particle separation 
vanishes. Unlike a truncated multipole expansion, the calculation presented 
here is exact. Following the systematic development of the procedure for an 
arbitrary number of particles, we demonstrate it explicitly with the simplest 

non-trivial example, the three-body system. 
02.20.Sv, 03.65.Fd, 31.15.-p 
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I. INTRODUCTION 



A system consisting of charged particles gives rise to A^(A^ — l)/2 pairwise Coulomb 
interaction terms in its Hamiltonian. Since only the two-body problem (A^ = 2) can actually 
be solved exactly, conventional atomic physics methods view the complete system of particles 
at the outset as a conglomerate of independent two-body systems, adding the interactions 
between these independent particles in a second step. This approach amounts to selecting 
a suitable subset of the A^(A^ — l)/2 Coulomb terms for which a solution in terms of "si- 
multaneous two-body wave functions" can be given. For instance, the independent-particle 
model for atomic systems treats each electron as interacting primarily with the nucleus (or 
with the ionic core in the case of valence electrons). Each (valence) electron's position intro- 
duces an independent spherical direction. Just as in the familiar solution of the hydrogen 
atom, each electron contributes to the angular part of the total system's wave function a 
spherical harmonic Yim of the angles specifying its direction in space. In the next step, all 
the electron-electron interactions are calculated by expanding the corresponding separations 
into Legendre polynomials of the inter-electronic angles, thus yielding the familiar multipole 
expansion. In more general terms, the independent-particle model first selects a specific par- 
ticle #1 (on physical grounds, typically the nucleus or the ionic core) and solves for each of 
the remaining A^ — 1 particles the two-body problem {1, j}, j = 2, . . . , N, thereby providing 
a basis for expanding the total wave function. In the next step, the interaction between 
particles i and j is calculated by adapting the reciprocal of their separation, l/r^j, to the 
coordinate system pertaining to the "two-body" basis functions for i and j {i, j = 2, . . . , N). 

The Coulomb interaction is singular whenever an inter-particle separation vanishes. Since 
these singularities are isolated from one another, they do not pose fundamental difficulties 
in calculating the Hamiltonian matrix. However, they give rise to cusps (discontinuous 
derivatives) in the wave function thus slowing down the convergence of partial- wave 
expansions of the wave function in their vicinity 0,^ . The independent-particle wave func- 
tions can only account (through the corresponding s-wave components) for the cusps arising 
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from vanishing separations rij, j = 2, . . . , N. Cusps due to vanishing Vij with i ^ 1 are not 
reproduced. Possible remedies to these shortcomings include the explicit use of coordinates 
like the inter-electronic angle 6ij or even Vij besides and r^; various approximations for 
l/vij 1^; or replacing ri,rj by r< = min(rj, r^), r> = max(rj,rj) in the Hamiltonian All 
these approaches amount to adapting the interaction operator l/rij to a single coordinate 
system. 

The present investigation explores the alternative approach of adapting the basis func- 
tions to match the relevant inter-particle separations Vij by utilizing several coordinate sys- 
tems simultaneously. We calculate a specific interaction term ~ ^/^ij in the coordinate 
system best suited for this particular purpose. Thus, the particle separations r^j dictate the 
choice of coordinate systems, the wave functions being transformed between the relevant 
reference frames to evaluate the different terms of the Hamiltonian. The success of this 
approach hinges on our ability to perform the numerous transformations between reference 
frames with high efficiency. Sec. ^ describes hyperspherical Jacobi coordinates appropriate 
for this task. Sec. |T| provides the main results by implementing the relevant transforma- 
tions for a system with an arbitrary number of particles and by constructing basis functions 
(harmonics) suitable for extensive transformation between reference frames. The resulting 
set of hyperspherical harmonics, derived here in the context of calculating the Hamiltonian 
matrix for a system of charged particles, has in fact much wider applicability. 

Hyperspherical coordinates and corresponding hyperspherical harmonics have been ap- 
plied in various areas of physics since the 1950s, for instance, in three-body scattering 
nuclear 0Jl^ and atomic [1TI|-[T9[] physics, as well as in quantum chemistry [^,0 . However, 



the sets of functions introduced in the present investigation are equivalent to, but much 
more flexible than, the hyperspherical harmonics discussed in P-pl[] . Beyond constituting a 
complete orthogonal set of functions appropriate for expansions, their frame-independence 
affords greater flexibility in analyzing selection rules and other relations between harmonics. 
These aspects are conveniently investigated through ladder operators; they are determined 
entirely by the symmetry properties of the A^-particle system, independent of any coordinate 



representation. Only calculating the non-vanishing matrix elements requires an appropri- 
ate coordinate representation of the generic basis functions. Sec. fV| illustrates the relevant 
procedures with the simplest non-trivial example, the three-body problem. The conclud- 
ing Sec. |V| discusses advantages and limitations of the new technique, as compared to the 
conventional multipole expansion. 



II. HYPERSPHERICAL JACOBI COORDINATES 

The inverse proportionality between pairwise interaction and inter-particle separation 
suggests replacing individual particles' positions with relative (Jacobi) coordinates. Starting 
from one pair of particles, the construction of Jacobi coordinates proceeds hierarchically by 
joining (the centers of mass of) increasingly complex groups of particles. This hierarchical 
structure is commonly referred to as a "Jacobi tree" |T0[. Alternative choices for the initial 
particle pair, as well as the ordering of successive particle groups, correspond to different 
Jacobi trees. In the next step, mass-scaling of the Jacobi coordinates fully exposes the 
symmetry of the kinetic energy operator for the multi-particle system. Let p, q denote two 
(groups of) particles with masses Mp, Mg and center of mass positions and Vg, respectively. 
Appropriate mass-scaling of the relative coordinate in the form ^ = {MpMg/ {Mp + 

-^g)}^^^!^?"''"?) removes the individual mass dependence from the expression for the kinetic 
energy: 

- y — K = — y ArPM- 2.1 

Setting the origin at the center of mass of the whole system allows discarding the CM's posi- 
tion and motion. With the individual mass factors removed, the kinetic energy (generalized 
Laplacian) displays complete symmetry under rotations in (3A^ — 3)-dimensional space. 

Hyperspherical coordinates exploit this symmetry by separating the "shape" of the sys- 
tem (described by 3A^ — 4 angular coordinates) from its overall "size" = J2k ^1 (with 
the dimension of a moment of inertia). Eigenfunctions of the Laplacian's angular part, 
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termed "hyperspherical harmonics" by generalizing the two- and three-dimensional cases, 
constitute a basis for expanding the complete wave function. In higher-dimensional spaces, 
the angular Laplacian's eigenvalues A(A + 3A^ — 5) are highly degenerate. In addition to 
the "grand angular momentum" A |p, the hyperspherical harmonics thus require numer- 
ous labels — analogues of the single "magnetic" quantum number m for three-dimensional 
spherical harmonics — for their unique specification. Different sets of hyperspherical harmon- 
ics, resulting from solving the Laplacian eigenvalue problem by separation of variables in 
alternative sets of coordinates, have been investigated extensively (see, e.g., 0,0,^] for 



systematic studies). However, their construction through separation of variables inevitably 
ties these harmonics to a specific coordinate system. They are thus not well suited for exten- 
sive transformation between different reference frames. The definition of harmonics based 
only on their behavior under the relevant transformations is the main result of this paper. 



to be derived in Sec. 111. 



The construction of each Jacobi tree starts with a pair of particles. Thus, each Jacobi 
tree contains at least one Jacobi vector joining two particles only (rather than centers of 
mass of particle groups). We denote such a vector as a "primary" coordinate. In any one 
Jacobi tree, up to \_N/2\ Jacobi vectors are directly proportional to actual inter-particle 
separations, Vij, the remaining relative coordinates necessarily involving larger complexes of 
particles. (Here and in the following, \_x\ denotes the largest integer not exceeding x.) We 
will calculate each of the A^(A^ — l)/2 Coulomb terms of the Hamiltonian in a Jacobi tree 
where it occurs as a "primary" coordinate. Obviously, all these terms then take the same 
form as the simple one-electron integral over 1/r in hydrogen, but the variables now refer 
to a multitude of different reference frames. Having thus eliminated all nested integrations, 
our next task consists in determining the transformations between different Jacobi trees. 

The most general transformation from one Jacobi tree to another (for the same system of 
particles) resolves into a sequence of elementary operations. Each elementary step consists 
in "transplanting" a sub-complex q from some particle complex {pq} to a complex {qr} 
it is achieved by a two-dimensional kinematic rotation p|] through the angle 



4> = tan~^ ,jMq{Mp + Mq + Mr)/ MpMr (2.2) 
in the (^p,^, ^p5,r)-plane of the (A^ — l)-dimensional space of mass-scaled Jacobi vectors: 

^'g r = COS (p ^p,g - sin ^pg^r (2.3a) 

$'p gr = sin (f) ip,q + cos $,pq^r- (2.3b) 

As each Jacobi vector is a vector in three-dimensional physical space, the basic rotation 
in the (^j, ^A;)-plane implies three rotations through the same angle in the (xj,Xfc)-, the 
{yj,yk)-, and the (zj, 2;jt)-planes, where {x,y,z) denote the Cartesian components of An 
elementary kinematic rotation through a finite angle in the (^j, ^fc)-plane then reads ^2 



Tij,iki'P) = exp(i0J^^.^J exp(i0J^^.^J exp(i0J^^.^J, (2.4) 

with infinitesimal rotation operators 
f d d\ 

= ( M— - t;— 1 , {u,v) = {xj,Xk), {yj^yk)A^j^Zk). (2.5) 

(For a heuristic explanation of ( p.4| ), recall the Taylor expansion of f{x + Ax) which may 
formally be written as exp{iAx p) f{x) with the infinitesimal "translation operator" p = 
—i-§^ corresponding to the quantum mechanical linear momentum (with ^ = 1). In the 
present context, the translation from x to x + Ax is replaced by a rotation through an angle 
from some direction in multi-dimensional space to a new direction.) In order to determine 
T^^.^^j,'s effect on the basis functions, we need to (i) express the generic rotation operators 

in terms of "ladder operators" whose action on the basis functions is straightforward to 
calculate, and (ii) construct basis functions suitable for such transformations. 

This section concludes with the outline of a strategy for determining the most efficient se- 
quence of two-dimensional rotations to achieve the transformation between arbitrary Jacobi 
trees. Note first that the individual terms in the Hamiltonian are referred to by means of 
particle labels, implying the need for appropriate antisymmetrization of the total wave func- 
tion under interchange of identical particles. This consideration dominates the construction 
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of the first Jacobi tree to wliicli all subsequent transformations are applied. In an atomic 
or molecular system, antisymmetrization concerns primarily the electronic part of the wave 
function. Thus the basic Jacobi tree starts out with a pair of electrons. It grows by adding 
one electron at a time, finishing with the addition of nuclei or ionic core(s). This procedure 



whose parentheses separate sub-complexes. Alternatively, one could start by first forming 
as many pairs of particles as possible, i.e., [N/2\ pairs, before joining these pairs into larger 
complexes. Although the antisymmetrization of an electron pair is particularly compact 
(namely, their relative angular momentum and their spin must add to an even value), anti- 
symmetrization among different pairs requires breaking up all these pairs, in essence going 
back to a canonical tree. Since breaking up sub-complexes involves about as many elemen- 
tary rotations as forming the new complexes, this procedure is not efficient. Nevertheless, a 
set of (for odd particle number) or A^ — 1 (even particle number) trees with the maximum 
number of [A^/2j pairs suffices to isolate all the A^(A^ — l)/2 inter-particle separations. A 
successful strategy therefore aims at building these particular Jacobi trees from the canonical 
tree. The simplest of these "pair-trees", (■ • ■ ((12)(34))(56) . . .), is obtained from the canon- 
ical tree with only [A^/2j — 1 rotations. In general, the pair {jk) with j < k is {k — j — 6j^i) 
rotations away from the canonical tree, but many other "useful" pairs are formed in the 
course of these "transplantations" . 



This section provides the tools required to rotate harmonics from one coordinate sys- 
tem to another. Our transformation matrices are analogous to the so-called Raynal-Revai 
coefficients The reader familiar with Lie algebra will of course recognize the relevant 
aspects of so{3N — 3), the algebra of rotations in (3A^ — 3) dimensional space. However, 



results in a "canonical" Jacobi tree []TU[, represented by the sequence of particle labels 




(2.6) 



III. TRANSFORMING THE BASIS FUNCTIONS 
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such familiarity is not presumed, and we hardly use the terminology of group theory. Our 
presentation reflects a more pragmatic point of view, emphasizing both the technical im- 
plementation as well as its relation to the physical application at hand rather than full 
mathematical generality. For the latter aspects, the reader should turn to the mathematical 
literature [Q. 



A. Rotations in d-dimensional space 

The most intuitive description of a rotation in three dimensions requires two elements: 
the (invariant) axis, and the angle of rotation. The rotation itself occurs in a plane perpen- 
dicular to the axis of rotation. In three dimensions, specifying the direction of the axis of 
rotation is equivalent to, but more economical than, describing the actual plane of rotation. 
For a higher dimensional space this is no longer the case, because there are several invariant 
directions perpendicular to a given plane. Thus, in d- dimensional space, a basic rotation is 
more appropriately characterized as occurring in the plane spanned by two coordinates rather 
than by an invariant axis orthogonal to it. The generic infinitesimal rotation operators then 
take the form (|2.5| ). The number of different planes, ^d{d — 1) in d dimensions, equals the 
number of basic rotations. Because rotations occurring in non-intersecting planes affect dif- 
ferent pairs of coordinates, they are evidently independent of each other; the corresponding 
rotation operators commute with one another. As there are \_d/2\ non- intersecting planes in 
a (i-dimensional space, the largest set of simultaneously commuting operators Juv contains 
i = ld/2\ elements. These commuting operators are commonly denoted as Hj, j = 1, . . . ,i 
2^ . For our application, the power of the Lie algebraic method derives from the structure 



and the properties of rotations being completely independent of any particular coordinate 
representation of the operators Hj. In fact, the following manipulations are most efficiently 
carried out in generic Cartesian coordinates without any reference to a specific Jacobi tree. 

Simultaneous eigenfunctions of the Hj with (integer) eigenvalues nij, j = 1, . . . , i, provide 
a basis for building hyperspherical harmonics suitable for extensive rotations between refer- 
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ence frames. If the infinitesimal rotation corresponding to Hj occurs in the (wj, Wj) -plane, 
the simultaneous eigenfunction of the £ Hj with eigenvalues rrij, respectively, reads 

F{{uk,Vk})yiXl{u,+iVjr^ (3.1) 
i=i 

with a function F that bears further specification in Subsect. [ill (J| . At this point, we only 
require that HjF = for all j. The set of eigenvalues /x = {mi,m2, . . . ,mi} serves as a 
label identifying different harmonics with the same value of A(A + d — 2) in the angular 
Laplacian's eigenvalue equation. Additional labels required for a unique specification will 
be introduced in Subsect. |III Q . 



Appropriate linear combinations of the remaining infinitesimal rotation operators act as 
ladder operators Ea satisfying 

[H,,E^]=a,E^, J = !,...,£, (3.2) 

with an ^-dimensional vector index ol with components Uj = 0, ±1. Here aj = indicates 
that Ea does not change the part of the eigenfunction pertaining to Hj, whereas aj = ±1 
means that it maps this part to the eigenfunction with eigenvalue ruj ± 1. The ladder 
operators interrelate eigenf unctions with different /x but degenerate A. Raising and lowering 
operators form Hermitian conjugate pairs E^^ = Ej^. 

In general (for i > 1), each ladder operator affects two of the rUj simultaneously; for 
odd-dimensional spaces a subset of i pairs of raising and lowering operators change one of 
the rUj only [^. Abbreviating the vector label o: by giving only its two non- vanishing 
components aj and ak, the ladder operator that raises ruj and simultaneously lowers ruk 
takes the form 

whereas the operator raising both rrij and rrik is represented by 



(Straightforward application of these operators to eigenfunctions (p.l|) verifies their behaving 

as ladder operators: they contribute to, or remove from, ( |3.1| ) factors {uj±ivj) and (uk^ivk) 

as appropriate for the intended "ladder operator action".) For odd-dimensional spaces, a 

residual coordinate, denoted here by w, does not occur in any of the Hj. The ladder operators 

changing only rrij (rather than a pair mj,mk) read then 

/ Q d \ 

Ef = -i {uj ± i Vj)- w— — — — — , j = 1, . . . , £. (3.3c) 

\ ow o[uj =F tVj) J 

(For d = 3, setting (u, w) = {x, y, z) and transforming the derivatives reveals the familiar 
pair of ladder operators ± ily, the single Hj occurring in this case coinciding with 1^.) 
This symbolic representation allows for efficient implementation on the computer. 

A subset of £ ladder operators E^. (and their Hermitian conjugates) suffices to interrelate 
all harmonics with the same eigenvalue A. A convenient choice p4| for these E^. is given 
by the ladder operators that raise mj and simultaneously lower mj+i for j = 1, 1. 
Ee/s form depends on whether d is even or odd. For d even, E^^ raises both m£_i and 
me, and for odd d, it raises me only, without changing any of the other mj. The set {sj}, 
j = 1, . . . ,i, then provides a basis for the ^-dimensional space of the vector labels ck and ^^. 
In an ^-component vector notation, this basis reads (we assume £ > 5 in order to expose the 
generic structure clearly) 

ei = (1, -1,0,0,. ..,0) 

62 = (0,1,-1,0,...,0) 

63 = (0,0,1,-1,. ..,0) 

: (3.4) 

e,_i = (0,0,..., 0,1,-1) 

(0,0,...,0, 1, 1) ioTd = 2i 
(0,0,..., 0,0,1) forrf = 2£+l. 
In general these basis vectors are not orthogonal in ^-dimensional space. Note, however, the 
following special cases: (i) In four-dimensional space (with i = 2), the two basis vectors ei = 
(1, —1) and 62 = (1, 1) are orthogonal, indicating that the two basic ladder operator pairs, 
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E±e^ and E±e2, commute with each other. In group theoretical language, this feature reflects 
the direct product structure 5*0(4) = 5*0(3) ® 50(3). However, the two 50 (3)-components 
do not refer to rrii and 1712 directly, but rather to (mi — 7112) /2 and (mi + 1712) /2. (ii) For 
rotations in three-dimensional space, we have ^ = 1, thus only one quantum number m which 
is being changed by one pair of ladder operators E^^^ = 1^. (iii) \n d = 2 dimensions, there 
are no ladder operators. Since all rotations occur in the same plane (the "only" plane of 
two-dimensional space), the order in which rotations through different angles are performed 
does not matter; they are all independent of each other. In the present formulation, the 
(single) rotation operator Hi generates all rotations. Each rotation is associated with its 
own harmonic function, exp(i0), the phase functions for different rotations (i.e., for different 
rotation angles cj)) not being related to one another through linear operators. 



B. Rotation of hyperspherical harmonics 

We now turn to the analysis of the transformation described by ( |2.4[ ). Note flrst that 
the three rotations occur in three non-intersecting planes, affording the more suitable rep- 
resentation exp(i0[Ja,^,,,^. -|- Jj^^.j^j.]) exp(i0 J^^.^^). Arranging the Cartesian components of the 
— 1 mass-scaled Jacobi vectors in the form 

{xi, yi, X2, 1/2, • • • , XN-1, Vn-i, Zi, ^2, • • • , zn-i}, (3.5) 

we choose the Hj by selecting pairs of coordinates from this list, starting from the left: 

Hj = Jx,y„ j = l,...,N-l (3.6a) 



Hn-i+u — Jz2k-iZ2k^ k — 1, . . . , 



N -1 



(3.6b) 



in the notation of (|2.5|) . We have stressed repeatedly our aim of deflning hyperspherical 
harmonics entirely through their behavior when acted upon by the operators Hj and E^., 
because these operators embody the structure and symmetry of the whole system in a frame- 
independent way. The explicit coordinate representation ( |3.6a| - PTFE| ) of the relevant opera- 
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tors appears at variance with our intended frame- independence. However, the spatial coor- 
dinates, {xj, yj, Zj), j = 1, . . . , N — 1, in ( p.6a| - |3l6ED should be viewed as generic Cartesian 
coordinates] they are not tied to any particular Jacobi tree. Arranging the generic Cartesian 
coordinates as we did in ( |3.5| ), on the other hand, does reflect physical considerations be- 
yond the purely mathematical structure of rotations in higher-dimensional space: The latter 
would instead label the coordinates most appropriately as (Xi, X_i, X2, X_2, . . . , Xi, X^i), 
supplemented possibly with an Xq for odd- dimensional spaces. Our arrangement takes into 
account that the dimension (3X— 3) arises from a product structure of an (X—1) -dimensional 
particle space with the three-dimensional physical space of each Jacobi vector. In particu- 
lar, the arrangement (|3.5| ) affords attributing relevant physical meaning to the first {N — 1) 
eigenvalues nij, (j = 1, ... ,X — 1): They represent the 2;-projections of physical angular 
momenta; their sum constitutes the z-projection of the total orbital angular momentum, 
an invariant of the system. 

As is evident from ( p.6aD , the coordinates {xj, Xk, yj, yk) occur in Hj and H^. The x- and 



y-parts of the rotation ( ^.41 ) thus involve the ladder operators that change mj and m^, only. 
A more detailed analysis shows that ( ^.3a| -c) can be inverted to read: 

JxjXk + JyjVk = ^jk + Ejk^ = Ea{jk) + E_a{jk), (3.7) 

for l<j<k<N — 1, with cx{jk) = J2^Zj ^s, i-e., the ^-component vector with +1 as its 
j-th component and —1 as its k-th component, all other entries being 0. 

The following consideration is central to our development, providing the crucial link 



between the ladder operator representation ( p.7|) for {Jx^x^ + JyjVk) their transformation 
matrix elements, by means of Euler-angle rotations. In analogy to the relation l^ = | (/+ + /_) 
familiar from rotations in three dimensions, we view the sum of a ladder operator Ea(jk) and 
its inverse (Hermitian conjugate) E_a{jk) as describing a rotation about an analogue of the 
X-axis . In three dimensions, an arbitrary rotation conventionally resolves into a sequence 
of three rotations: about the z-axis, about the resulting y'-axis, and about the new z'-axis, 
through the Euler angles {ip,6,ilj), respectively ||22|. In terms of these three Euler-angle 
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rotations, a rotation about the x-axis through an angle results from the following sequence 
of operations: The first Euler-angle rotation about the z-axis through the angle = — 7r/2 
rotates the y-axis onto the original x-axis; in the second step one rotates about this y'-axis 
(which is the original x-axis) through the angle 6 = (p; the third Euler-angle rotation finally 
moves the x'-axis back to the original x-direction by rotating about the z'-axis (lying in the 
original (t/2;)-plane) through the angle ip = 7c/2. Wigner's (i-symbol is the matrix element 
for the rotation of a spherical harmonic about the y-axis, d^^'mi^) — (^m'| exp(i0/j,)|Y/m), 
whereas the two 2;-type rotations only contribute phase factors exp(2mV/2 — imTi/2) p2| , 
giving for the rotation of a spherical harmonic about the x-axis 

eM^<PQ\yim) = e''"'^ ^11^(0) e"^"^^- (3.8) 

Generalization of the concept of Euler-angle rotations from three to more dimensions 
hinges on the following key observations: m and m' serve to distinguish between degenerate 
harmonics with the same I. In higher-dimensional spaces, the vector indices and /x' play 
the same role. However, the parameter I in the d-symbol indicates not only the angular 
Laplacian's eigenvalue Z(/ + 1), but also — more importantly — the range of possible m values, 
— / < (m, m') < I. More precisely, it sets the upper limit of 21 for the number of times either 
one of the two ladder operators /+, can act in direct succession before necessarily mapping 
any spherical harmonic Yjm to zero. (Since l+Yim ^ Yi,m+i, it takes I — m < 21 "raising 
operator" steps to reach Yu from Yim- Applying Z+ in succession {21 + 1) times necessarily 
raises m beyond its upper limit. The absence of a corresponding harmonic translates to 
{l+y''~^^Yim = 0.). Viewed in this way, the harmonic's parameter m indicates its "position" 
along the string (of length 2/ + 1) of degenerate harmonics interrelated by a ladder operator. 
We now extend these concepts to rotations in more than three dimensions. 

The single m (Z^'s eigenvalue) in three dimensions belongs to a string ranging from — / to 
+1, accessed by the single pair of ladder operators In more than three dimensions, we 

replace it with an ^-component vector fj, whose components are changed (typically in pairs) 
along different strings labeled by corresponding ladder operator pairs Ea and E^a- (Recall, 
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however, that the ^ vectors provide a basis for the vectors a; thus any ladder operator 
Ea can be expressed as products of E^..) The role of m as an indicator of the harmonic's 
position along the single string in three dimensions extends therefore to higher dimensions 
if we project the "indicator" onto the "direction" of any particular string ex. of interest, 
i.e., if we define 



rria = , = , (3.9) 

a OL oc - a 

where the scalar product in the denominators accounts for ck's two non-zero components 
in the case of a pairwise change of m^'s. In analogy to / in the three-dimensional case, 
a parameter Aq, determines how many times the ladder operators E±ct can be applied in 
direct succession. The modulus of nia in (|3.9| ) provides a lower limit for the relevant "string 
length" Act- Since the ej form a basis for the ex., we anticipate that Ac, emerges from an i- 
component vector A = (Ag^, . . . , A^J. The latter provides the additional parameters required 
for a unique specification of harmonics (besides /j, and the eigenvalue A) . Note that in the 
general case Aq, is not identical with the eigenvalue A, at variance with the situation in three 
dimensions. 

To summarize the procedure so far, we note that the x- and y-parts of the finite rotation 
(2.4) affect only the two components mj and rrik of the vector fj,, shifting rrij by an integer 



amount n while simultaneously changing rrik by the same amount in the opposite direction. 
The projections ( p.9|) evaluate to 

rriaUk) = ^{rrij - rrik)] (3.10a) 



m 



OL{jk) 



-{m'j - m[) = mo,{jk) + n, (3.10b) 



with non-vanishing matrix elements for the rotation of a hyperspherical harmonic ^a./x.a 
into another harmonic Yp^i ^j^i x occurring only if and /x' lie along the same a-string of 
harmonics with degenerate A, i.e., if an integral number n of ladder operator steps E±a(jk) 
separates /x' from /j,. With our previous considerations on Euler-angle rotations leading to 
( p.8|) , this matrix element reads 
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where we omitted a's parameters j, k for brevity. The Kronecker symbols 6 ensure that only 
harmonics with the same eigenvalue A, lying along the same a-string, are connected. The 
factor 2 multiplying the rotation angle reflects the absence of the factor 1/2 in the RHS 
of ( |3.7| ), as compared to the expression for Ix = ■^{1+ + thereby effectively doubling the 
rotation angle in (^.8|). 

It remains to apply the same concepts to the z-part of the rotation {T^). This part 
affects only ^t's components that correspond to the involving Zj and z^. Because we 
chose to gather all the ^-coordinates after all (x?/)-pairs in (|375| ), the Jacobi vector indices 
j, k are shifted relative to ^t's components for the ^-coordinates. To simplify the notation, 
we use 



i = N-l + 



J + 1 



N -1 + 



k + 1 



(3.12) 



so that and now refer to the operator and quantum number pertaining to the coordi- 
nate Zj, with the same connection between H^, m^, and Zk- Depending on the indices j and 
/c, i.e., depending on the positions of the coordinates Zj and Zk in the sequence ( |3.5| ), Jz^z^ 
takes alternative forms in terms of ladder operators or if^- The possible cases are: 



1. k even and j = k — 1. According to ( |3.6bD , Jz^zk coincides with -ffAr_i+fc/2 = H^^. The 
rotation matrix element is simply 



(>A',M',A'| exp(i0 J^^.^J|Ya,m,a) = exp( 



'^A',A'^A',A'5/x',/x- 



(3.13) 



2. k even and j < k — 1 even. We find 

Jz,z, = I {e:- + i?r/ + + e;~) . (3.14) 

Viewing again the sum of a ladder operator and its inverse (Hermitian conjugate) 
as proportional to an analogue of a rotation about the x-axis, we readily recognize 



15 



the last expression as the sum of two x-type rotations (compare with ( |3.7| )'s central 
expression). With the explicit formulae ( p.3a| - [3l3BD , it is straightforward to verify 



that the ladder operator E'^~ raising and lowering commutes with the "raising- 
raising" operator -Ej^^. The two x-type rotations are therefore independent of each 
other. As in the previous discussion of the x- and y-parts, an integer number of steps 
E'^~ must connect (m^, m^) to {m[, m'^), i.e., m[ = mi^ + ni, m'^ = — rii. However, 
the second rotation involves the same components of /x, only this time changing both 
in the same direction: m[ = + n2, m'l^ = + n2- Compatibility of the two 
conditions requires rii = n2 = 0, thus 

= ^("^t - ^k) = m'f3_(,^) (3.15a) 
= ^{m, + m«) = m;3^(,^), (3.15b) 

and the extra phase factors occurring in (p.ll| ) drop out in this case: 

(^A',/x',A'| exp(i0 JzjZk)\YA,fJ.,x) = dmlz]mp_ (0/2)(im^+,''m0_^ (0/2) (5a',a5a',A (3-16) 

The factor 1/2 multiplying the rotation angle stems of course from the factor 1/4 in 



3. k even and j < k — 1 odd. In this case. 



J.^,, = i - + - e;,-) , (3.17) 



obviously the sum of two y-type rotations. The "selection rules" for and are 
the same as in the previous case, and since the extra phase factors distinguishing an 
X- from a y-type rotation cancel in ( p.l6| ), the matrix element is identical to the one 
in case 2. 

4. k < N — 1 odd and j even. JzjZk turns into 

JzjZk = — [EttT - - + E~~) . (3.18) 



Ai 
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Exactly the same considerations apply again, except that the difference of two y-type 
rotations corresponds to a rotation through the angle —0/2 for the Ef^ part. 

5. A; < — 1 odd and j even, yielding 

Jzr^u = \ {El' + - Et,^ - E;-) , (3.19) 

i.e., the difference of two x-type rotations. The matrix element coincides with the one 
of case 4. 

6. k = N — 1 odd and j even, does not occur in any of the H^, and the relevant ladder 
operators act on only: 

= I {Et + e;) . (3.20) 

This is an x-type rotation, exactly as in three dimensions. Projecting /j, and fi' onto 
the appropriate direction 7(i): = m^, m'^^^-^ = m[ = m~^i^o + we obtain for the 

matrix element 

(1a',m'.a'I exp(z0 J.,.J|1a,m.a) = ^^!^^,L^(0)e'("'^-'"-)"/'5A',A5A',AV,M+--y (3-21) 

7. A; = — 1 odd and j odd. 



i {e+ - E;) , (3.22) 



a y-type rotation. Its matrix element differs from the one in case 6 only by the absence 
of the phase factor exp(iz/7r/2). 



Thus, the 2;-part of the rotation ( p.4|) takes essentially three different forms: (i) a phase 
factor exp(zmK0), diagonal in the /x, if Jz^z^ coincides with if^; (ii) an x- or y-type rotation 
through the angle 0, changing into m[ = + z/, as in three dimensions, if Zk is the 
unpaired coordinate zn-i of an odd-dimensional space; (iii) the product of two rotations 
through the angles 0/2 and (—1)^^0/2, respectively, but diagonal in ^t, if Zj and z^ belong 
to different if^ and H,^. 
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C. Harmonics suitable for rotation 



In the previous subsection, we have completely determined the matrix elements for the 
rotations of our interest when applied to hyperspherical harmonics, before actually specifying 
these functions explicitly. This was possible because we expressed the matrix elements in 
terms of the labels ( "quantum numbers" ) (A, /x. A) identifying the harmonics, rather than 
through integrals in coordinate space. So far, the harmonics are functions of the generic 
Cartesian coordinates {xj,yj,Zj). Upon rotation from one Jacobi tree to another, these 
coordinates transform into another set {x'j , y'j , Zj) , j = 1, . . . , N — 1, but the harmonics, 
when expressed in terms of the new coordinates, retain their functional form. In this sense, 
these hyperspherical harmonics are frame-independent. 

Direct solution of the angular Laplacian's eigenvalue problem by separation of the 
(3A^ — 4) coordinates in the second-order differential equation leads to hyperspherical har- 
monics represented by standard spherical harmonics and Jacobi or Gegenbauer polynomials, 
depending on the choice of hyperspherical coordinates |1T0|JT7| , |2T|| . However, these harmonics 
are not simultaneous eigenfunctions of all the Hj. Consequently, they are not suitable for 
our purpose, because rotating the harmonics requires knowledge of the ladder operators' 
effects on these functions, at least for the base set of ladder operators E±^^, j = 
which in turn requires uniquely specifying the harmonics with labels (A,^i,A). We now 
construct complete sets of functions defined exclusively by their behavior under the action 
of the first-order differential operators Hj and -Eiej, J = 1? • • • , ^- The resulting functions are 
"hyperspherical harmonics" , too, because they satisfy the generalized Laplacian's symmetry 
under rotations. 

The hyperspherical description separates the "hyperradial" momentum from the gener- 
alized angular momentum. Each of the Cartesian coordinates Xj,yj, Zj, {j = 1, . . . ,N — 1), 
is proportional to the "hyperradius" R. The angular Laplacian's eigenvalue parameter A 
determines the harmonic's degree by setting the radial scale as R^. Neither the Hj nor the 
ladder operators E±(:^ affect this radial factor, as is to be expected of angular momentum-like 
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operators. With (3A^ — 4) angular coordinates, a complete specification of the harmonics 
requires (3A^ — 5) labels in addition to A. The vector /x provides £ = [(3A^ — 3)/2j of them 
in the form of the eigenvalues rrij of all the Hj. The remaining labels are taken from the 
^-component vector A consisting of the "string lengths" Ag^ for the ladder operators E±ej. 
The first {£ — 2) components of A suffice to reach a total of (3A^ — 4) labels if is odd, as 
do the first (£ — 1) components for even. 

It remains to determine the function F introduced in (|3.1|) . The requirement Hj F = for 
all j implies that F depends only on {uj+Vj) (and possibly on the single unpaired coordinate 
w in case (3A^ — 3) is odd). Furthermore, if J2j = A the product Y[j{uj + i Vj)"^^ already 
accounts for the radial factor R^, i.e., F = const, in this case. This circumstance suggests 
constructing complete sets of degenerate harmonics with the same A as follows: We set mi 
to its maximum value and all other rrij = 0, i.e., /x = (A, 0, . . . , 0) [Q. For j = 1, . . . ,i, the 
modulus of the projection m^^ = /x ■ ej/^ej ■ sj) sets the lower limit for the corresponding 
string length Ag^. However, Ag^ cannot be larger than m^^ either, for applying any of the 
raising operators -E+g^ necessarily results in a vector fi' = fi + Sj having J2j l^^'jl > i-G-; 
in a harmonic not belonging to the same set of degenerate functions. With A completely 
specified, the first harmonic reads 

Ya,^,x = CA (xi + t y,)\ A = (mg„ 0, . . . , 0); ^ = (A, 0, . . . , 0), (3.23) 

where ca denotes a normalization constant. Starting from this first function, we generate 
harmonics by applying all the lowering operators E_^^ first to (|3.23|) , then applying the E_^. 



to the harmonics so obtained, and so on in a recursive procedure: 

>^A,M-..A' = ^ ^-,>A,.,A (3.24) 

'(Ag, +mgJ(Ag, -mg, + 1) 



for all Gj-strings that have not yet terminated, i.e., the strings having m^^ > — Ag^.. The 
recursive nature of the process suggests gathering together the harmonics in "levels" , with 
the level of a harmonic indicating the number of lowering operator steps separating it from 
the first harmonic (|3.23|) . At level with the harmonic (|3.23|) , the ei-string is the only string 
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with non-zero length, and we can only generate one harmonic of level 1. But for i > 1, this 
level- 1 harmonic will already have several non- vanishing m'^, and corresponding A'^^ , giving 
rise to additional e^-strings starting from this level. Repeating this procedure level by level, 
we work our way down until finally reaching the last harmonic with /x = (—A, 0, . . . , 0). The 
procedure stops automatically, because applying any of the lowering operators to this last 
harmonic maps it to zero (exactly as IJYi-i = in three dimensions). 

In d dimensions, the total number of different hyperspherical harmonics with the same 
"grand angular momentum" A is [^, p. 265] 



, „ 2A + d-2/A + d-2\ 

For a given A, the procedure outlined above generates exactly dim(A, d) independent func- 
tions, i.e., a complete set of hyperspherical harmonics. However, some of these harmonics 
still require modifying for our purpose, as the following observations illustrate. 

Note first that with our choice for E^., the first {^ — 1) ladder operators always change 
two TTij in opposite directions, thus leaving I'^jl = A invariant when starting from the 
first harmonic. For all these harmonics, the function F reduces to a constant. At some 
stage during the above construction, however, an e^-string will appear along which the sum 

l^jl no longer remains constant. Since the ladder operators do not change the overall 
radial factor i?^, the function F accounts for any "missing powers" in Wjiuj + ivj)"^^ for 
harmonics having fi with J2j l^il < A. For instance, in three dimensions the lowering 
operator /_ {E^ in ( |3.3c| )'s notation) removes one power of (x + iy) from the "highest" 
harmonic r'-Ya ~ (x + iyY = r' sin' Oexp^ilcf)) while simultaneously adding a factor F = z = 
rcosO to the next harmonic r^Yi^i^i ~ 2; (x + iyY~^ = r' cos^sin'"""^ ^exp(i[/ — In the 

same way, repeated application of /_ generates the higher-degree Legendre polynomials in 
cos 6 making up F in this case. 

So far, A and the vector /x provide enough information to uniquely specify the hyper- 
spherical harmonics; we need their additional label A only for the purpose of determining 
transformation matrix elements, not to distinguish the harmonics from one another. For 
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these harmonics, the construction described above provides the appropriate function F. 

In a higher-dimensional setting, however, the E^^ do not all commute with each other. 
It is thus possible to arrive at different functions with the same along different ladder 
operator sequences, starting from the first harmonic ( |3.23| ). This is exactly the situation of 
"degenerate eigenvalues": in this case, /x, the set of eigenvalues of the Hj, has multiplicity 
higher than one. Although our construction generates the appropriate number of indepen- 
dent functions, thus providing a basis system for the higher-dimensional eigenspace associ- 
ated with /i, we now need additional labels — to be taken from A — to distinguish the different 
eigenfunctions with degenerate eigenvalue /x. Suppose, therefore, the vector ^t occurs with 
multiplicity z/ > 1. Our procedure generates v independent functions $p, p = 1, . . . , i/, all 
having the same Y{j{uj + ivj)"^^ . These functions differ only in their F. Harmonics suitable 
for the calculation of rotation matrix elements are expressed as linear combinations 

V 

= E«-P'^P' (3-26) 
p=i 

and the requirement 

= 0, for J = 1, . . . , w (3.27) 

determines the coefficients a^p for a = 1, . . . ,1/ (and thus the functions F). In this way, ap- 
propriate sets of parameters (A^^ , • • • , ^e^^^^ ) serve to distinguish the harmonics by enforcing 
specific lengths for the different ej-strings. As noted previously, the number of additional 
labels required is jmax = ^ — 2 (or £ — 1) for even (or odd)-dimensional spaces, respectively. 
With this modification, even the harmonics corresponding to degenerate vectors /j, show the 
desired behavior under rotations in d dimensions, and their transformation matrix elements 
can be deduced from their "quantum numbers" (A, ^t, A) directly. 

Finally, a remark concerning the degeneracy of ^t seems in order. In three dimensions, 
it is impossible to arrive at the same = m along different strings, because there is only 
one pair of ladder operators. Nevertheless, even in this case we note that the m-components 
arising for a given / (= A) occur again for /' > /. Due to the mostly pairwise change of 
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TUj in higher-dimensional spaces, we expect that a vector fi occurring for a given A will 
appear again as a label for harmonics with A' = A + 2, A + 4, . . .. For d > A there are 
non-commuting ladder operators. Because the number of different pathways leading from 
the highest harmonic to a specific /x increases with the length of this path, /x's multiplicity 
increases with A. Furthermore, since different functions F imply contributions from different 
(«! + f|)-terms, higher multiplicity — even in spaces with more than three dimensions — can 
only arise for A — J2j l^j l > 2. A detailed analysis of the recursion in A with /x held fixed 
confirms these expectations, yielding for the multiplicity of n 

muh(^; A, d) = K + ^l (3.28) 



with p 





d-3 


; Q = 


2 



2 

a useful result to test the implementation of the (recursive) procedure. 

IV. EXAMPLE: TWO-ELECTRON SYSTEM 

For the purpose of illustration, we apply the method outlined in the preceding sections 
to a Coulombic three-body system. Specifically, we consider a two-electron atom or ion 
with atomic number Z, i.e., a system with only one heavy particle. The case of diatomic 
molecules (two-center Coulomb system) requires additional modifications of the hyperspher- 
ical approach ^7 . 



A. Coulomb interactions among three particles 

After elimination of the CM motion, a three-particle system requires two Jacobi vectors 
for its description. We use three different Jacobi trees, ^j, r]i, <^j, {i = 1,2). With Vi, r2, and 

denoting the positions of the two electrons and the nucleus, respectively, the relevant 
mass-weighted Jacobi vectors read 



-Yiri-r,), (4.1a) 
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m = y j^J^^M ~ "^i) = ^^^-^ ^1 + si^/^ ^2; (4.2) 



Ci = - ^n) = COS7 6 + sm7 6; (4.3) 



where 



-^--V 2(m4) ' ^^-^^^ 

. / M;v + 2Me 

""^ = -V 2(M. + M.) ' 

and 7 = — /3. M^r and Me denote the nuclear and electron mass, respectively. The vectors 
T72 and will not be needed, since the Coulomb interaction among the three particles takes 
the form 



1 M^M, Z Z \ 

At this point, the familiar approach using the transformation of the coordinates would ex- 
ploit ([4.2| - [4.3|) in a multipole expansion of the electron-nucleus interactions in terms of the 



coordinates ^1 and ^2- In our method, however, we apply the (kinematic) rotations (|42- 
[4.3| ) to the wave functions instead. This amounts to calculating the three interaction matrix 
elements in three different coordinate systems. The immediately obvious advantage of this 
method is that each of the three terms takes exactly the same form. Each of the integrals 
reduces to the same radial integral as in the textbook example of hydrogen. Higher-order 
multipoles and nested integrations over powers of the coordinates do not occur. 

B. Choice of Hj and ladder operators 

The kinetic energy has rotational symmetry in six dimensions. Partitioning the six coor- 
dinates into pairs defines three non-intersecting planes, and thus i = 3 mutually commuting 
rotations. We choose the corresponding first-order operators Hj as Jxiyi, Jx2y2i Jz^z2- 
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Here {xi,yi, zi,X2,y2, ^2} denote the Cartesian components of the mass-weighted Jacobi 
vectors. Accordingly, the harmonics — simultaneous eigenfunctions of the three Hj with 
eigenvalues n = (mi, m2, m^) — take the form 

2 

F{xl + yf, xl + yl zf + zj) Hixj + tsign{mj) y^)'™^' (^1 + ^sign(m3) ^2) (4.6) 

i=i 

where we have rewritten ( p.l| ) so as to yield non-negative powers of the hyperradius for 
either sign of the rrij. The appropriate base set of three ladder operator pairs E±^^ is then 
specified by the vector labels ei = (1, —1, 0), £2 = (0, 1, —1), and £3 = (0, 1, 1), corresponding 
to ^2^3^, and ^21^ in ( CT^ ). 



C. Labeling the basis functions 

Five angles specify each point on a sphere of fixed hyperradius R in six-dimensional space. 
Consequently, our hyper spherical harmonics require five labels — related to the numbers of 
nodes in the five angular variables — for their identification. Besides A and /x = (mi, m2, m-^) 
we need one more label, A^^ , the string length along the ladder operator sequences spanned 
by E±ei- 

The set of harmonics so defined differs from the more familiar set labeled by quantum 
numbers {Ii,mi,l2,m2,na) [|I7|,^. However, the latter set of harmonics does not make use 



of the third eigenvalue m3. Since rotating the harmonics is accomplished by acting on them 
with ladder operators which in turn act on all the rrij, the conventional harmonics specified 
by {li,mi,l2,^2,na) are not suited for our purpose. They may, of course, still serve as a 
basis set in an application, being then expanded into the new set labeled by (A, /x, A^J prior 
to the actual rotation. A straightforward transformation links the two basis sets. 



D. Explicit construction of harmonics with degenerate /x 



As an example, we derive the expressions for some harmonics that are not completely 
characterized by A and /x. Consider for instance harmonics with A = 4, /x = (1,0,1). 
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According to ( |3.28|) , there are two harmonics with these labels, to be distinguished by 



the additional parameter A^^. With the particular vector fi of this example, we obtain 



m 



|, thus setting the lower limit for Ag^. The two essentially different ways of arriving 



at the set of labels (1, 0, 1) from harmonics of the next-lower level are given by an E^i-step 
from fi' = (2,-1,1), and by an E^^,_^-step from fi" = (1,1,2). Both of these labels have 
multiplicity 1 since J2j \ = A; the function F{. . .) in the corresponding harmonics reduces 
to a normalization constant. Up to these normalizing factors, Ci and C2, ( [4.6|) gives these 
harmonics as 

^4,(2 -1,1), f = ci {xi + iyif {x2 - iy2){zi + ^^2), (4.7a) 
^4,(1,1,2), 1 = C2 (xi + iyi) (x2 + iy2){zi + iz2f. (4.7b) 

Applying the appropriate ladder operators to these harmonics provides a basis for the two- 
dimensional eigenspace of degenerate /x = (1, 0, 1): 

$1 = E_g,F4^(2,-l,l),f 

= C3 (xi + iyi){zi + iz2) {[xl + yl] -2[xl + yl]) (4.8a) 

$2 = E_e3y4,{l,l,2),l 

= C4 (xi + iyi){zi + %Z2) {2[xl + yl] - [zl + zl]), (4.8b) 

with constant factors C3, C4 to be ultimately absorbed into the normalization. Because $1 is 
obtained from (2 _i 1) 3 by applying -E_g^, it behaves under the relevant rotations exactly 
as is required for ^4 (1 1) Apart from the normalization constant C3, we thus find 

^4,{i,o,i),f = C3 {xi + iy{){z^ + iz2) {[xl + yl] - 2[x^ + yl]). (4.9) 

However, the second harmonic with the same /x = (1, 0, 1) does not simply coincide with 
$2, because the latter has the same string-length Ag^ = | as $1: One easily verifies that 

Eg,$2^0. (4.10) 

We can, however, construct a harmonic (1 1) ^ ^ linear combination of $1 and $2 by 
requiring 

25 



E,,{a^i + b^2)=0, (4.11) 
yielding the condition Sacs — 26c4 = and thus 

^"4,(1,0,1)4 = C5(a:i + iyi)izi + iZ2) {2[xl + yl] + 2[xl + yl] - ?,[zi + zl]). (4.12) 
The two harmonics with "degenerate" ^ = (1, 0, 1) are now distinguished by their respective 

E. Rotating the harmonics 

Eq. ( [4.2|) describes the coordinate transformation between the Jacobi trees (^1,^2) and 
{'Hii'n^)- Accordingly, the hyperspherical harmonics transform as 



|A,M,A,,)^ = ^d;^-)(/5) |A,/x',Ae,)^ (4.13) 



where the subscripts on the ket vectors indicate the respective Jacobi tree. The transforma- 
tion matrix elements are given by 

^mV^I/^) = (A,^', AeJexp(i/5[J:,i^2 + Jyiy2 + -^^i^2])|A,M, Aei) 

= (A, A,, I expitPlE,, + E.,, + H^])\K, fx, A,,), (4.14) 



according to Subsec. [1IIB| . Following the analysis given there, we find 



D^^^l\f3) = e^[-'-l-/2rfL^i(2/3) 5^, ^^^^^^ (415) 

with 

m = 2^^'^ ~ "^2) (4.16a) 

m' = ^{m[ - m'^) (4.16b) 

n = m[ — mi = —m'2 + m2 (4.16c) 

in terms of ^-components. Replacing the angle (3 with 7 = —j3 yields the expressions for 
the transformation to Jacobi tree 
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F. Symmetry properties of basis functions 



Due to the very high degeneracy of harmonics with the same grand angular momentum 
A as expressed in ( p.25| ), it is important to exploit various symmetries of the functions in 
order to reduce the size of the hyperspherical basis. Reflection through the origin of the 
coordinate system transforms all six coordinates into their negatives. According to ( [4.6| ), 
the harmonic |A, /x, Ae^) picks up a factor (— i^l^il+l^al+lmai ^j^^^der this operation. Since all 
the ladder operators change two of the rrij at a time, the sum in the exponent has the same 
parity as A. This first observation thus restricts the basis set by allowing only even A for 
even-parity states and odd A for odd-parity states. 

Consider next the harmonics' symmetry under interchange of identical particles (i.e., the 
two electrons). This interchange is most easily described in the Jacobi tree {^1,^2}: 

^ 6 ^ 6- (4.17) 

Accordingly, the harmonic |A, mi, 1712, ms, Ae^) turns into the harmonic |A, mi, m2, — ms, Ag^) 
under interchange of the electrons while acquiring a factor (— i)™i+™3+s' (^-^^here S denotes 
the total spin of the two electrons). Antisymmetrized basis functions may thus be labeled 
by non-negative only, and = is possible only for harmonics with (mi + S) even. 

Furthermore, when choosing the first-order operators Hj we have arranged the six co- 
ordinates in a way that allows us to identify mi and m2 with the 2;-projections of the 
three-dimensional relative angular momenta Zi and Z2. Therefore, the sum mi + m2 = mtot 
is the z-component of the coupled (total) orbital angular momentum L. The fact that the 
Coulomb interaction does not couple states with different mtot reduces the size of a basis 
consisting of antisymmetrized harmonics accordingly. This point is particularly interesting 
because the system's invariant mtot restricts our basis, even though the corresponding total 
orbital angular momentum L is not defined in this basis. The reason for this seemingly 
surprising fact lies in our use of first-order operators only; hence Lz can be identified, but 
not the second order operators Z| or L^. The absence of the invariant three-dimensional L 
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is the main trade-off we have to accept when treating all transformations as rotations in a 
genuinely six-dimensional space especially Sec. 10.2., p. 267ff|. 



G. Calculation of matrix elements 



While Cartesian coordinates prove most appropriate for manipulating the hyperspherical 
harmonics using ladder operators, hyperspherical coordinates lend themselves for the calcu- 
lation of matrix elements. Specifically, the familiar representation of Cartesian coordinates 



xi = -Rcosasin6'i cos(/9i, 
Ui = -Rcos a sin ^1 sin 
Zi = -R cos a cos 6'i , 



X2 = R sin a sin O2 cos Lp2 , 
1/2 = R sin a sin 62 sin ip2 , 
Z2 = -R sin a cos 6^2, 



R^F sinl^^l eie'"^"^' sinl'^^l ^se*™^^^ (cos a cos 61 + isign(m3) sin a cos 62 



(4.18) 



transforms (14. 61) into 



\l™.3| 



(4.19) 



where F is a function of (cos^ a sin^ 6*1), (sin^ a sin^ 62), and (cos^ a cos^ 9i + sin^ a cos^ 62). 
Note that this form remains the same, regardless of whether the Cartesian components 
(x, y, z) refer to the Jacobi vectors in tree ^, in 77, or in The relevant interaction operator 
entering into the matrix element is always l/(/2cos a), for each of the three pairwise Coulomb 
interactions, with the angle a referring to a different coordinate system in each case. The 9- 
and a-integrals arising in the calculation are related to the Euler Beta function [EH] , namely, 



2-K 



de sin^e cos'' 9 



:i + (-i: 



(p-l)!!(g-l)!! 





7r/2 



{p + q)\\ 



-c, 



pq 



da sin^ a cos'' a 







(p-l)!!(g-l )!! 
(p + g)!! 



(4.20a) 
(4.20b) 

(4.20c) 



where 



p, q both even 



otherwise. 



(4.21) 
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A multipole expansion would have forced us to split the last integral into two parts, with 
different integrands depending on whether | cos jS^i \ is greater or smaller than | sin (3^2 1 • The 
resulting integral could only be expressed as a sum of factorial quotients, rather than a single 



term, as in ( |4.20c| ) 



Gathering together all the pieces developed in this section, we find for the angular part 
of the matrix elements 

(A', A',jFc|A, A,,) = (A', m', K,\- ^ - — |A, /X, A,,) (4.22) 



X 



X (A', /xi, a;J^|A, /X2, a,,), (4.23) 
^ cos a 

with the further simplification 7 = —(3. The double summation (transforming the bra and 
ket vectors between Jacobi trees) seems to spoil the present approach's advantage over a 
multipole expansion. After all, the latter also leads to two summations, namely, a sum 
over the multipole order and another summation stemming from the analytical evaluation 
of the nested integral over a. Note, however, that the present method achieves significantly 
more with two summations: it accounts for all the cusps in the wave functions whenever an 
inter-particle separation vanishes. 

Finally, the particular case of a three-body system involves two relevant Jacobi vectors 
only. All possible Jacobi trees are thus related to one another by rotations in the same plane 
(^1, ^2)- This particularity of the three-body problem might suggest arranging the Cartesian 
coordinates of the Jacobi vectors in the following way: 

{xi,a;2,?/i, 2/2, ^1,^2}, (4.24) 

rather than our choice ( p.5| ). The Hj resulting from the above arrangement of coordinates 
coincide with JxiX2: Jyiy2i ^"^^ JziZ2, thereby simplifying the rotation of harmonics: All 



rotations reduce to the first case of (ziZ2)-type rotations discussed in Sec. [1IIB| , with matrix 
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elements 

(A', /x', A',J exp(z/5[J,,,., + Jy,y, + J.,.J)|A, /x, A,,) = c^^'^^+'^^+'^^^^'aVm^a^.a,, • (4.25) 

This orthogonality relation virtually eliminates the double summation over /xi, /X2 in ( [4.23^ . 
However, the corresponding harmonics are now functions of (xi ± 1x2)-, {yi ± W2), and 
{zi ± 22:2). Transforming to hyperspherical coordinates using 

xi = -Rsinai cosa2 cos(y9i, X2 = -Rsinai cosa2 sin(y9i, 

Hi = i?sinai sina2 COSV92, 1/2 = -R sin a 1 sin 0:2 sin ^92, (4.26) 

zi = -Rcosai cosy^s, Z2 = -Rcosai siny^s, 

yields Hj = —i-^, with < < 2n and < < n/2. While the Hj, as well as 
the harmonics, now obviously attain their simplest form, the choice ([4.24|) has two serious 
drawbacks: (i) None of the three eigenvalues mj of the Hj have physical significance, and (ii) 
Vl^il = (^1 + 2/1 + ^i)'^^"^ is a (complicated) function of all five angles (ai, 0^2, v^i, V52, V's)- 
The latter problem is solved by using hyperspherical coordinates ( [4.18| ) instead of ( [4. 261) , 
yielding for each of the combinations (xi ± 1x2), {yi ± ^2)5 and {zi ± 2^2) a sum of two 
terms (as opposed to the single terms obtained for {x,j ± iyj) in the previous sub-sections). 
Expanding powers of these binomials leads to two additional summations, leaving us with 
no net gain. 

V. CONCLUSION 

Recognizing the independent-particle model's failure to account for cusps due to variables 
on which the wave function does not depend explicitly, we have developed a method that 
satisfies Kato's cusp condition through reference frame transformations. By transforming 
the wave function to the appropriate reference frame, we expose the wave function's cusp 
arising from vanishing of any given particle separation rij. In addition to satisfying the cusp 
condition on the wave function, this technique also simplifies the calculation of the pairwise 
Coulomb interaction ~ l/rj^, as compared to the conventional multipole expansion. 
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Implementation of the approach outlined above resolves into three major tasks, all ad- 
dressed in the present investigation: (i) the systematic study of the relevant transformations 
between reference frames; (ii) the definition of functions suitable for such transformations; 
(iii) the determination of the transformation matrix. While the literature on Lie algebra 
provides ready-made solutions to problem (i), it usually fails to do so for (ii) and (iii). Fur- 
thermore, the mathematical literature does not exploit the particularity of an atomic or 
molecular A^-body system. 

More specifically, problem (i) is solved by using mass-scaled Jacobi coordinates, since the 
transformations between reference frames reduce then to generalized rotations in (3A^ — 3) 
dimensions. The Lie algebra so{3N — 3) describes these transformations completely, em- 
bodied in the sets of commuting rotation operators {Hj} and ladder operators {E±^.}, 
j — 1, . . . , L(3iV — 3)/2j. To solve problem (ii) mentioned above, we have introduced basis 
functions defined entirely through their behavior under infinitesimal rotations, i.e., when 
acted upon by the operators Hj and E±^.. Simultaneous eigenfunctions of all the Hj con- 
stitute appropriate basis functions, classified further according to their matrix elements for 
the transformation between reference frames. The latter step removes possible ambigui- 
ties whenever the Hj have degenerate eigenvalues. Finally, by extending the concept of 
Euler-angle rotations from three to higher dimensions, we have provided the solution to task 
(iii), the determination of the transformation matrix elements. For an arbitrary A'^-body 
system, each change of reference frames considered here resolves into a sequence of basic 
transformations described by these matrix elements. 

In an A'"-body system, the high dimensionality arises from the product structure of the 
(A^ — 1) -dimensional particle-space and the three-dimensional physical space. The basic step 
in the transformation between reference frames reduces to a two-dimensional rotation in a 
plane of particle-space. Upon expansion of the particle-space variables into their physical- 
space components, the basic rotation induces three two-dimensional rotations in (3A^ — 3)- 
dimensional space that arc analyzed using the Lie algebra so{3N — 3). To take advantage of 
the efficiency offered by the Lie algebraic method, we treat the transformations as rotations 
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in a genuine (3A^ — 3)-dimensional space. Thus, we work with the first- order diSerential op- 
erators Hj and E±^. throughout. Arranging the generic Cartesian components {xj,yj,Zj) of 
the particle-space variables appropriately, we can partially recover the product structure 
characteristic for the physical apphcation at hand: The eigenvalues oi {N — 1) among the 
Hj carry physical significance and we interpret them as ^-projections of orbital angular mo- 
menta. Their sum represents the total angular momentum's ;s-projection, an invariant of the 
system. However, neither the individual angular momenta, nor the coupled (total) angular 
momentum can appear in our treatment, as they are represented by second-order operators. 
Other methods exploit the product structure of the (3A'" — 3)-dimensional coordinate space to 
a larger extent by solving the Laplacian's (second-order) eigenvalue problem through separa- 
tion of variables. However, the resulting hyperspherical harmonics incorporate these features 
of the three-dimensional physical space at considerable cost: Transforming these functions 
between reference frames proves very inefficient, requiring essentially their expansion into 
the equivalent sets of harmonics introduced in the present investigation. 

In conclusion, we have demonstrated by explicit construction the possibility of describing 
a system of N charged particles without recourse to multipole expansions. Kato's cusp 
condition is satisfied implicitly through reference frame transformations of the wave function. 
We introduced appropriate basis functions and discussed their symmetries under particle 
interchange and reflection through the origin. We derived the matrix elements for the 
transformation between reference frames for arbitrary numbers of particles, and we showed 
that the interaction matrix elements attain a simpler form as compared to the conventional 
multipole expansion. Because the formalism underlying our approach relies exclusively on 
first-order differential operators, it does not incorporate angular momenta, thereby limiting 
its usefulness for bound state problems. However, the technique introduced here provides 
significant simplifications in scattering problems where partial angular momenta are not 
resolved. 
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